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It is shown that the concept of a fireball with a definite filamentary structure naturally 
emerges from the analysis of the spectra of Gamma- Ray Bursts (GRBs). These results, 
made possible by the recently obtained analytic expressions of the equitemporal surfaces 
in the GRB afterglow, depend crucially on the single parameter Ti, describing the effective 
area of the fireball emitting the X- and gamma ray radiation. The X- and gamma ray 
components of the afterglow radiation are shown to have a thermal spectrum in the 
co-moving frame of the fireball and originate from a stable shock front described self- 
consistently by the Rankine-Hugoniot equations. Precise predictions are presented on a 
correlations between spectral changes and intensity variations in the prompt radiation 
verifiable, e.g., by the Swift and future missions. The highly variable optical and radio 
emission depends instead on the parameters of the surrounding medium. The GRB 
991216 is used as a prototype for this model. 

Keywords: black hole physics — gamma rays: bursts — gamma rays: theory — ISM: 
clouds — ISM: structure — radiation mechanisms: thermal 



1. Introduction 

It is generally agreed that the GRB afterglow originates from an ultrarelativis- 
tic shell of baryons with an initial Lorentz factor 70 ~ 200-300 with respect to 
the interstellar medium (see e.g. Ruffini et alP, Bianco & Ruffim^ and references 
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therein). Using GRB991216 as a prototype, in Ruffini et alM^we have shown how 
from the time varying bolometric intensity of the afterglow it is possible to infer 
the average density {riism) — 1 particle/cm'^ of the IntcrStellar Medium (ISM) in a 
region of approximately 10^^ cm surrounding the black hole originating the GRB 
phenomenon. It was shown in Ruffini et al."^ that the theoretical interpretation of 
the intensity variations in the prompt phase in the afterglow implies the presence 
in the ISM of inhomogeneities of typical scale 10^^ cm. Such inhomogeneities were 
there represented for simplicity as spherically symmetric over-dense regions with 
("iirn) — (uism) Separated by under-dense regions with (nff^) ~ 10~^ {nism) 
also of typical scale ^ 10"'^^ cm in order to keep {riism.) constant. The goal of the 
present letter is to show how the GRB spectral observations allows to reconstruct 
the full 3-dimensional ISM ffiamentary structure with (nll''^ ~ 'R.~^ ("-iiL)' where 
the factor TZ describes the effective area of the fireball emitting the X and gamma 
radiation. 

The main results are: 

a) The X and gamma radiation in the afterglow has a sharp and stable structure 
with evolution toward the hard and the soft part of the spectrum describable by 
the Rankine-Hugoniot equations. Such radiation has a thermal spectrum in the co- 
moving frame of the fireball and originate from the shock front created ahead of 
the ultrarelativistic shell of baryons generating the afterglow. 

b) Necessary condition for such a process to occur is a filamentary structure in the 
ISM in order to explain the observed luminosity and to fulfill the necessary opacity 
condition. 

c) The knowledge of the intensity variation of the prompt radiation (Ruffini et 
al.El) , of the luminosity in fixed energy bands in the afterglow (Ruffini et al.l^ and 
of the variability of the spectra (presented here) are then sufficient to determine 
self-consistently the 3-dimensional distribution of the ISM filamentary structure. 
As in previous articles, we use GRB 991216 as a prototype for our model. 

2. The structure of the shock front 

The interaction between the ultrarelativistic shell of baryons and the ISM into 
which it is expanding is mediated by a shock front with three well-defined layers 
(see e.g. sees. 85-89, 135 of Landau & LifsitJ3, ch. 2 and sec. 13-15 of Zel'dovich 
& Rayzei® and sec. IV, 11-13 of SedoJ^. From the back end to the leading edge 
of this shock front there is: 

a) A compressed high-temperature layer, of thickness A', in front of the relativistic 
baryonic shell, created by the accumulated material swept up in the ISM. 

b) A thin shock front, with a jump AT in the temperature which in the comoving 
frame can be described by the Rankinc-Hugoniot adiabatic equations: 



AT ~ (3/16)mp(5uVfc - 1-5 x lO" [5v / (IQ^kms^^)] ^ K , 



(1) 
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where 6v is the velocity jump, is the proton mass and k is Bohzmann's constant, 
c) A pre-shock layer of ISM swept-up matter at much lower density and tempera- 
ture, both of which change abruptly at the thin shock front behind it. 

At larger distances ahead of the expanding fireball the ISM is at still smaller 
densities. The upper limit to the temperature jump at the thin shock front, given 
in Eq.(^, is due to the transformation of kinetic energy to thermal energy, since 
the particle mean free path is assumed to be less then the thickness of the layer 
(a). The thermal emission of the observed X- and gamma ray radiation, which as 
seen from the observations reveals a high level of stability, is emitted in the above 
region (a) due to the sharp temperature gradient at the thin shock front described 
in the above region (b). 

The optical and radio emission comes from the extended region (c). The de- 
scription of such a region, unlike the sharp and well defined temperature gradient 
occurring in region (b) , requires magnetohydrodynamic simulations of the evolution 
of the electron energy distribution of the synchrotron emission. Such analysis has 
been performed using 3-D Eulerian MHD codes for the particle acceleration mod- 
els to produce the energy spectrum of cosmic ra ys a t supernovae env elop e fronts 
(see e^ McKee & Cowid^, Tenorio-Tagle et allHI, Stone & Normaiff^l^ Jun & 
JonellSI), Other challenges are the magnetic field and the instabilities. We mention 
two key phenomena: first, the importance of the development of Kelvin- Helmholtz 
and Rayleigh- Taylor instabilities ahead of the thin shock front. The second is the 
dual effect that the shock front has on the ISM initial magnetic field, first through 
the compression of the swept-up matter containing the field and secondly the am- 
plification of the radial magnetic field component due to the Rayleigh- Taylor insta- 
bility. Simulations of both effects (see e.g. Jun 

& JoneJEI 

and references therein), 

modeling the synchrotron radio emission for an expanding supernova shell at vari- 
ous initial magnetic field and ISM parameters, shows for example that the presence 
of an initial tangential magnetic field component may essentially affect the result- 
ing magnetic field configuration and hence the outgoing radio flux and spectrum. 
Among the additional effects to be taken into account are the initial inhomogeneity 
of ISM and the contribution of magnetohydrodynamic turbulence. 

In the following we focus uniquely on the X- and gamma ray radiation, which 
appears to be conceptually much simpler than the optical and radio emission. It 
is perfectly predictable by a set of constitutive equations, which leads to directly 
verifiable and very stable features in the spectral distribution of the observed GRB 
afterglows. In line with the observations of GRB 991216 and other GRB sources, 
we assume in the following that the X- and gamma ray luminosity represents ap- 
proximately 90% of the energy fiux of the afterglow, while the optical and radio 
emission represents only the remaining 10%. 
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Distance from the black hole (cm) 

Fig. 1. The ratio Mb /Mism is plotted as a function of the radial coordinate of the expanding 
baryonic shell. 



3. The theoretically predicted spectra of GRB 991216 

For GRB 991216 the initial conditions at the beginning of the afterglow era are: 
To = 1.9428 X lOi"* cm, 70 ~ 310.13, the baryonic matter Mb ^ 1.6122 x 10^° g, 
Etot ^ 4.8298 X 10^3 erg (Ruffini et alP). By fitting the relative intensity of the 
Proper GRB (P-GRB, see Ruffini et aM) and of the afterglow (Ruffini et alP), 
the average density has been shown to be slowly varying: (riism) — 1 particle/cm^ 
{{nism) — 3 particles/cm"^) in the early (late) afterglow phases (Ruffini et alPj). 
The integration of the Taub relativistic Rankine-Hugoniot equations of motion for 
the baryonic matter has been given in Bianco & RuffinPI The MbJ-S the "piston" 
in the language of Landau & LigitJ3, Zel'dovich & Rayze Sedov^. To ensure the 
validity of the treatment, it is important to check that the ratio Mb /Mism, where 
Misrn = (4/3)7rmpnisT?i('"'^ — f^), is always much bigger than 1 in the entire range 
of the observed afterglow, which is indeed the case for r < 4 x 10^^ cm (see Fig.0. 

The observed intensity variations in the "prompt" emission have been shown to 
originate in the over-dense and under-dense ISM anisotropics in Ruffini et al.El. 

Here the fundamental new assumption is adopted (see also Rufhni et al.EI) that 
the X- and gamma ray radiation during the entire afterglow phase has a thermal 
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Fig. 2. The temperature in the co-moving frame of the shock front corresponding to the den- 
sity distribution with the six spikes A,B,C,D,E,F presented in RufRni et al.^. The dashed line 
corresponds to an homogeneous distribution with Uigrn = 1- 



spectrum in the co-moving frame. The temperature is then given by: 
T, - [A^;i„t/(47rr2ATa7^)]'/* , 



(2) 



where Ai?i,it is the internal energy developed in the collision with the ISM in a time 
interval At in the co-moving frame, a is the Stefan-Boltzmann constant and 



n^A,ff/A, 



(3) 



is the ratio between the "effective emitting area" of the afterglow and the surface 
area of radius r. In GRB 991216 such a factor is observed to be decreasing during 
the afterglow between: 3.01 x 10"^ > 7^ > 5.01 x lO^^^ (RufRni et eM). 

The temperature in the comoving frame corresponding to the density distribu- 
tion described in Ruffini et al."^ is shown in Fig. |21 

We are now ready to evaluate the source luminosity in a given energy band. The 
source luminosity at a detector arrival time f^, per u nit solid angle dil and in the 
energy band [vi, 1^2] is given by (see Ruffini et al.ll^: 



dE. 



7 



dtd;dn 



EQTS 



^WCOSI^A ^ -^W {vi,V2,Tarr)dY., 



(4) 
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where Ae — AEint/V is the energy density released in the interaction of the ABM 
pulse with the ISM inhoniogeneities measured in the comoving frame, A = 7(1 — 
{v I c) cos ■&) is the Doppler factor, W (yx , , Tarr) is an "effective weight" required to 
evaluate only the contributions in the energy band \y\ , 1^2] , is the surface element 
of the EQTS at detector arrival time on which the integration is performed (see 
also Ruffini et alP) and Tarr is the observed temperature of the radiation emitted 
from dS: 

Tarr = T,l [7 (1 - (vlc)cOS'd) (1 + z)\ . (5) 

The "effective weight" W {i'i,i'2, Tarr) is given by the ratio of the integral over 
the given energy band of a Planckian distribution at a temperature Tarr to the 
total integral aTarr- 



W {vi,U2,Tarr) = 




where p {Tarr, ^) is the Planckian distribution at temperature Tarr'- 

p {Tarr, 1^) - [2 / h^) hv / [e^- / i'^^.^^) - l) (7) 

In Ruffini et alAwe applied this theoretical framework for the analysis of the ob- 
served luminosities in fixed energy bands (50-300 keV and 2-10 keV) of GRB991216 
over an observation time lasting 10^ s and we gave, as well, a self-consistent ex- 
planation of the observed hard-to-soft transition of the GRB spectra and of their 
variability during the GRB afterglow. We here illustrate the expected spectral dis- 
tributions for selected values of the photon arrival time at the detector, during the 
early phases of afterglow (see Fig. . These predictions appear to be of the great- 
est interest in view of the forthcoming Swift^ mission and future high-precision 
missions. The slopes of the four spectra appear to be constant with a = 0.9 in 
the range 1-5 keV, while they appear to be extremely varying in the energy range 
40-80 keV. We have, in fact, over an arrival time interval of a few seconds (see Fig. 
EJ): /3 = -0.4 at = 16.26 s, /? = -0.7 at if = 18.98 s, /3 = -0.25 at if = 22.73 s. 
It is particularly attractive to analyze observationally the very precise correlation 
predicted in our model between the intensity variations during the prompt emission 
and the effective power-law indexes of the high-energy and low-energy parts of the 
spectra. 

The most important feature is that, although the original emission in the co- 
moving frame is thermal, the final theoretically predicted spectra appear to be 
highly non-thermal! The final spectra, in fact, originate from a superposition of an 
infinite series of contributions along the EQTS depending on the angle between 
the emission region and the line of sight (see Ruffini et al.l'^. Each contribution 
carries with itself a distinct thermal spectrum boosted by a specific Lorentz factor. 
In order to visualize the departure from a thermal spectrum, we have also shown 
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Fig. 3. Theoretically predicted spectra (solid curves) for selected values of the photon arrival 

time at the detector: if = 3.488 s (a), = 16.26 s (b), = 18.98 s (c), = 22.73 s (d). The 
dotted straight lines corresponds to the slopes a in the low (/3 in the high) energy bands at 1-5 
keV (40-80 keV). It is interesting to compare and contrast the solid curves, representing the actual 
spectral distribution, to the corresponding boosted black body spectra at = (dashed curves, 
see text). 

in each figure the expected spectrum from an ideal emission process in which all 
the energy radiated is assumed to originate at = with the same Lorentz gamma 
factor, which is of course a perfect blue-shifted black body spectrum. 

4. On the physical origin of the 7?. factor and of the thermal 
radiation 

Wc arc now left to explain the physical origin a) of the thermalization process and 
b) of the TZ factor. Both these problems are solved if we assume that the ISM 
matter in the over-dense region is distributed in a filamentary structure in which 
each filament has an average density 

{n{L) = {<L)lT^. (8) 
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which is then varying from {n{ll^^ ~ 10^ {nism) at r ~ 1.0 x 10^® cm up to 

nlll^ — 10^'^ {nism) at r ~ 10^^ cm. Since the X- and gamma-ray emission only 
occurs in the filamentary structure, it follows from Eq.ljHl that the effective emitting 
area is indeed the one given by Eq. (PJ . 

The filamentary structure, which has the typical length of an over-dense region 
Ai? ~ 10^^ cm, is optically thick: 

T^aTin{:l)^R>\, (9) 

where ctt is the Thomson cross-section. The internal temperature Tint of the com- 
pressed layer of thickness A' in the above region a), measured in its co- moving 
frame, is given by: 

r„t = {AiSint/ [l&^r^i^ATon (A'/Ai?)] , (10) 
Since we must have Tint > 7s, we then have: 

A' < Ai?/ (47^) . (11) 

5. Conclusions 

Having conceptually explained both the thermalization process and the TZ. factor 
as originating from the filamentary structure, the "inverse problem" is now open to 
the solution: to infer, from the observed spectra and intensity variation, the actual 
three dimensional extension of the filamentary structure. 

The existence of filamentary structures will indeed sharpen the description of 
the intensity variations presented in Ruffini et alP'^ This will in principle allow to 
describe all the spikes in the observed prompt radiation, as well as the event D 
of Ruffini et al.l^', which were not explainable in the average spherically symmetric 
approximation. 
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